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ABSTRACT 


The  main  motivation  for  this  thesis  study  is  that 
significant  workload  for  aging  transport  aircrafts  is 
related  to  dent  removal  from  fuselages.  This  thesis  is  a 
preliminary  investigation  of  aircraft  fuselage  dents  using 
the  Finite  Element  Method  (FEM)  via  FEA  ABAQUS  software.  We 
investigated  single  impact  dent  on  fuselage  panel  at  various 
locations  and  impact  speeds.  The  material  used  for  our 
finite  element  models  is  Aluminum  Alloy  2024-T3,  a  typical 
material  used  for  fuselages  of  older  transport  aircrafts. 
Our  finite  element  model  consisted  of  impact  analysis, 
buckling  prediction  analysis,  and  postbuckling  analysis 
successively.  These  analyses  were  performed  on  both 
stiffened  and  unstiffened  aluminum  panels.  We  found  that, 
depending  on  dent  status  in  aluminum  panel,  the  dent  may 
increase  or  decrease  buckling  load  of  the  panel  compared  to 
that  of  the  virgin  panel  (undented)  .  The  buckling  load  of 
panel  with  low  velocity  impact  is  generally  lower  than  that 
of  the  virgin  plate.  As  the  impact  velocity  is  increased, 
buckling  load  of  dented  panel  increases  exceeding  buckling 
load  of  virgin  plate.  In  addition,  we  also  noticed  an 
existence  of  critical  impact  velocity  at  which  the  buckling 
load  of  the  dented  panel  reached  maximum  and  after  which 
will  start  to  decrease. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

With  today' s  rapid  increase  in  number  of  aging 
aircraft,  the  aerospace  industry  is  fighting  to  reduce 
costs  on  maintenance  and  repair.  Life  extension  of 
aircraft  structures  has  become  a  major  focus  of  the 
industry.  The  two  main  costly  problems  which  have  caught 
the  attention  of  the  scientific  community  are  the  fight 
against  fatigue  and  corrosion  [Ref.  3]  .  Much  studies  have 
been  and  are  actively  being  conducted  on  crack  fatigue, 
multiple  site  damages  (MSD) ,  aircraft  skin  inspection 
technology  and  amongst  others  [Ref.  4]  .  All  these  studies 
are,  in  one  way  or  another,  related  to  fatigue  and 
corrosion  crack  damages  of  aircraft  fuselage.  Although  a 
tremendous  workload  is  also  being  spent  on  dent  repairs  and 
removal  in  large  transport  aircrafts  [Ref.  1],  not  much 
study  have  been  done  on  investigating  the  static  stability 
of  dent  or  impact  damages  on  fuselage  panels. 

After  an  extensive  literature  search  and  review,  only 
a  few  studies  have  been  found  to  involve  impact  damages  on 
airplane  structure.  One  in  particular,  by  Cornells  Guijt 
from  US  Air  Force  Academy,  performed  experimental  and  some 
computer  modeling  study  on  fuselage  dents.  The  study  used 
a  special  impact  swing  hammer  to  develop  various  dent  sizes 
on  A1-2024-T3  fuselage  panel.  Although  this  paper  focused 
mainly  on  fatigue  investigation,  it  did  present  some  static 
stability  data  for  comparison. 
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B.  SCOPE  OF  RESEARCH 

This  paper  is  a  preliminary  study  on  fuselage  dents. 
The  ultimate  goal  is  to  be  able  to  provide  guidelines  on 
fuselage  dent  repairs.  This  paper  studies  the  effect  of 
fuselage  dents  on  compressive  failure  loadings  using 
computer  simulation.  The  fuselage  dents  analyses  were 
modeled  using  a  computer  simulation  program  called 
ABAQUS/CAE  for  both  unstiffened  and  stiffened  panels.  After 
the  deformed  fuselages  were  created,  they  were  subjected  to 
compressive  loading  and  eigenvalue  calculation  analyses. 
The  buckling  load  results  between  the  compressive  failure 
analyses  and  eigenvalue  analyses  were  compared.  The 
buckling  behaviors  of  unstiffened  and  stiffened  panels  were 
also  presented. 
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II .  MODELING 


A.  UNSTIFFENED  FUSELAGE  PANEL 

The  unstiffened  fuselage  panel  was  modeled  using 
ABAQUS/CAE  as  a  0.508  m  by  0.508  m  (20  in.  by  20  in.)  3-D 
deformable  shell  panel.  The  shell  panel  was  created  via 
extrusion  of  one  of  the  x-axis  edge.  The  material  property 
was  defined  as  elastic-plastic  with  shell  thickness  of 
0.0015875  m  (0.625  in.)  .  Density  of  the  shell  plate  was 
2780  kg/m3,  elastic  modulus  was  73.1  GPa  with  Poisson  ratio 
of  0.33.  An  80x80  S4R  (quadrilateral  shell)  finite  element 
mesh  was  used  for  the  shell  plate.  A  summary  of  the 
material  property  and  geometric  dimension  are  tabulated  in 
Table  1&2  and  a  finite  element  mesh  of  the  fuselage  panel  is 
given  in  Figure  1 . 

Table  1.  Material  Property  of  Shell  Plate  A1-2024T3 

(Ref.  5) 


Young's  Modulus  of  Elasticity 

73.1 

GPA 

Yield  Strength 

3.45E+08 

Pa 

Ultimate  Strength 

4 . 83E+08 

Pa 

Mass  Density 

2780 

Kg/m3 

Elongation  at  Failure 

15% 

Shell  Plate  Dimension  (LxWxt) 

0.508x0.508x0.003175 

mxmxm 

Table  2.  Geometric  Dimensions  of  Fuselage  Panel  (after 

Ref.  1) 


Model  Geometry 

Unstiffened  Plate  Dimension 

0.508x0.508 

m 

Shell  Thickness 

0.15875 

cm 

Stiffened  Plate 

Dimension 

0.508x0.508 

m 

Shell  Thickness 

0.0015875 

cm 

Stiffener  Dimension 

1 . 905x1 . 905x1 . 905 

cm 

Shell  Thickness 

0.15875 

cm 
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Figure  1.  Finite  Element  Mesh  of  Shell  Panel- (element 

size=0 . 635  cm) . 

B.  STIFFENED  FUSELAGE  PANEL 

The  Stiffened  fuselage  panel  was  modeled  using 
ABAQUS/CAE  similar  to  the  unstiffened  panel.  The  main 
difference  was  the  addition  of  hat  shaped  stiffener 
supports.  The  0.508  m  by  0.508  m  (20  in.  by  20  in.) 
stiffened  panel  consisted  of  a  set  of  hat  stiffeners  at 
20.32  cm  (8  in.)  apart.  Each  of  the  joint  between 
stiffener  and  main  panel  was  modeled  to  be  twice  as  thick 
as  (0.3175  cm  or  0.1250  in.)  the  main  panel  itself.  The 
mesh  size  was  the  same  as  the  unstiffened  panel  but 
included  the  mesh  for  the  two  stiffeners.  The  height  and 
width  of  winged  portion  of  each  stiffener  were  both  1.905 
cm  (0.75  in.)  .  The  material  property  used  was  same  as 
listed  in  Table  1.  A  detail  dimension  of  stiffened  panel 
is  shown  in  Table  2. 
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Figure  2.  Finite  Element  Mesh  for  Stiffened  Panel 

(element  size=0.635  cm)  . 

C.  UNSTIFFENED  FUSELAGE  PANEL  WITH  PROTRUDING  DEFORMATION 

An  additional  study  was  conducted  to  investigate  the 
stiffening  effect  caused  by  the  present  of  fuselage  dents  on 
compressive  failure  loads.  The  model  was  created  using 
ABAQUS/CAE  similar  to  the  unstiffened  fuselage  panel.  The 
model  consisted  of  a  0.508  m  x  0.508  m  deformable  shell 
panel  with  protruding  hollow  cone-like  structure  as  shown  in 
Figure  3.  The  deformable  shell  panel  with  20.32  cm  (8  in.) 
diameter  center  cutout  was  created  separately  from  the 
protruding  structure.  ABAQUS  has  a  feature  in  which  parts 
with  same  uniform  material  property  can  be  merged  or 
combined  into  one  piece.  In  this  case,  the  cutout  shell 
panel  was  merged  with  the  protruding  structure  to  form  our 
final  modeling  part.  The  material  property  used  for  the 
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protruding  structure  was  same  as  the  rest  of  the  shell  panel 
(Table  1)  .  Figure  4  showed  a  side  view  of  a  protruding 
structure  with  the  dimensional  property. 


Figure  3 . 


Finite  Element  Mesh  of  Shell  Panel  with 
protruding  structure  with  20.32  cm  (8  in.) 
diameter  cutout. 
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Figure  4.  Side  view  of  protruding  structure. 
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III.  SIMULATION 


A.  INTRODUCTION  TO  ABAQUS/CAE 

All  computer  simulations  in  this  study  were  conducted 
using  a  finite  element  simulation  program  called  ABAQUS. 
ABAQUS  is  an  advanced  finite  element  analysis  that  provides 
complete  and  powerful  solutions  for  linear  and  nonlinear 
engineering  problems.  It  is  a  suite  of  finite  element 
analysis  modules  consisting  of  preprocessor,  solver,  and 
postprocessors.  Preprocessor  allows  the  creation  and 
assembly  of  modeling  parts.  Their  solver  package  includes 
ABAQUS /Explicit  and  ABAQUS/Standard .  ABAQUS  Viewer  is  the 
postprocessor  which  allows  visualization  and  post  processing 
of  analysis  results.  Their  most  complete  package, 
ABAQUS/CAE,  is  a  fully  interactive  environment  equipped  with 
graphical  user  interface  (GUI)  of  menus,  icons,  and  dialog 
boxes.  It  provides  the  complete  working  environment  for 
users  to  create  models,  submit  jobs  for  analysis  via  the 
solver  package,  and  view  and  process  results  [Ref.  2] . 

B.  MODELING  AND  SIMULATION  PROCEDURE 

Both  the  stiffened  and  unstiffened  plate  models  were 
subjected  to  three  different  analyses.  The  first  analysis 
was  the  dynamic  impact  simulation  in  which  an  artificial 
dent  is  created  via  a  5.08  cm  (2  in.)  diameter  rigid  ball  at 
various  locations  and  impact  velocity.  The  deformed  plates 
were  than  imported  along  with  the  material  states 
(stress/strains)  using  ABAQUS/CAE  in  order  to  undergo 
compressive  failure  analysis,  which  was  conducted  in  two 
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different  ways.  The  first  study  was  the  incremental 
compressive  loading  to  the  impacted  panels  until  the 
buckling  of  the  panels.  In  this  study,  the  residual  stress 
and  strain  caused  by  the  impact  was  included  in  the  model  as 
well  as  the  dent  shape.  The  other  study  was  undertaken  for 
eigenvalue  analysis  of  the  dented  panels.  In  the  latter 
case,  the  deformed  shapes  of  the  panels  were  considered  in 
the  analysis,  but  not  the  residual  stress  and  strain.  Both 
results  were  compared  to  the  buckling  loads  of  the  virgin 
panels  before  impact. 

A  total  of  21  impact  simulation  runs  were  conducted  and 
each  took  average  of  15  hours.  Each  of  the  28  compression 
test  simulations  took  from  1-3  hours  to  complete.  Finally, 
the  time  required  for  each  eigenvalue  calculation  analysis 
was  approximately  30  minutes.  All  analysis  were  conducted 
using  a  Dell  Dimension  Desktop  equipped  with  single  Pentium 
4  processor  at  2.4  GHz  and  1  GB  RAM.  A  summary  of  all 
simulations  is  given  in  Table  3. 

C.  DYNAMIC  CONTACT  IMPACT  SIMULATION 

The  ABAQUES/Explicit  analysis  engine  was  used  to 

perform  the  dynamic  impact  simulation  for  each  fuselage 

model.  ABAQUS/Explicit  possessed  a  powerful  nonlinear 

dynamic  analysis  capability  which  allowed  us  to  define 

contact  interaction.  In  this  analysis,  surface-surface 

contact  interaction  was  defined  for  the  rigid  ball  impactor 

and  the  deformable  fuselage  panel.  The  normal  behavior  of 

interaction  property  was  defined  as  hard  contact.  Contact 

pair  separation  after  impact  was  allowed  in  order  to 

eliminate  contact  load  after  impact  on  the  deformed  fuselage 

panel.  The  frictionless  formulation  was  used  for  the 
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tangential  contact  interaction  behavior.  The  Kinematics 
Contact  Method  was  used  as  the  mechanical  constraint 
formulation.  The  boundary  conditions  of  the  fuselage  panel 
were  rigidly  fixed  on  all  four  edges  with  an  initial 
condition  of  the  model  defined  for  the  rigid  ball  impactor 
velocity.  An  enhanced  hourglass  mesh  control  was  chosen  in 
order  to  reduce  the  effect  of  hourglass  stiffness.  This 
mesh  control  property  propagated  into  subsequent  analyses 
and  did  not  require  further  adjustment.  For  each  fuselage 
panel  model,  a  series  of  impactor  velocity  was  used  to 
investigate  the  various  dent  sizes.  Total  simulation  time 
step  was  set  at  1  second.  A  complete  detail  of  the  creation 
of  this  analysis  model  is  provided  in  Appendix  A. 

D.  BUCKLING  PREDICTION/EIGENVALUE  CALCULATIONS 

ABAQUS/Standard  analysis  engine  was  used  for  buckling 
prediction  analysis.  Individual  deformed  geometry  created 
from  the  impact  analyses  was  imported  into  a  new  model  to 
obtain  eigenvalues  for  the  first  5  modes.  A  simple  support 
boundary  condition  was  used  on  the  two  opposite  ends  while 
the  other  two  opposite  ends  remained  free.  The  in-plane 
compressive  loading  was  applied  to  one  of  the  simply 
supported  edges.  A  uniform  shell  edge  load  of  1  Newton  was 
used.  The  total  predicted  buckling  load  is  found  by 
multiplying  the  respective  eigenvalue  value  by  1  Newton. 
The  eigensolver  selected  for  this  analysis  was  Subspace 
iteration  method.  A  detail  description  of  the  creation  this 
model  is  provided  in  Appendix  A. 
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E.  COMPRESSION  FAILURE  ANALYSIS/POSTBUCKLING  SIMULATION 

ABAQUS/Standard  analysis  was  also  used  for  compression 
failure  analysis.  Unlike  the  eigenvalue  calculation 
analysis,  the  deformed  fuselage  geometry  along  with 
respective  residual  stresses  and  strains  were  imported  into 
a  new  model  in  order  to  perform  the  compression  test 
simulation.  Boundary  conditions  used  were  simple  supports 
as  described  in  the  previous  section.  An  incremental 
uniform  shell  edge  load  was  used  for  all  compression  test 
analyses.  Static  Riks  Iteration  Method  was  chosen  for 
buckling  calculation. 


Table  3.  Fuselage  Panel  Simulation  Summary. 


A1-2024-T3 

Fuselage  Panel  Simulations 

Panel  Type 

Impact  Type 

Impact 

Location 

Impact 

Speed 

(m/s) 

Eigenvalue 

Calculation 

Compression 

Test 

CJndeformed  Panel 

N/A 

N/A 

N/A 

Yes 

No 

Undeformed  Stiffened 

Panel 

N/A 

N/A 

N/A 

Yes 

No 

Deformed  Unstiffened 

Panel 

Ball 

(r=l " ) 

Center 

10 

Yes 

Yes 

Ball 

(r=l " ) 

Center 

30 

Yes 

Yes 

Ball 

(r=l " ) 

Center 

35 

Yes 

Yes 

Ball 

(r=l " ) 

Center 

55 

Yes 

Yes 

Ball 

( r-1 " ) 

Center 

60 

Yes 

Yes 

Ball 

(r=l " ) 

Center 

65 

Yes 

Yes 

Ball 

(r=l " ) 

Center 

70 

Yes 

Yes 

Deformed  Stiffened  Panel 

Ball 

(r=l " ) 

Center 

10 

Yes 

Yes 

Ball 

(r=l " ) 

Center 

30 

Yes 

Yes 

Ball 

(r=l " ) 

Center 

35 

Yes 

Yes 

Ball 

(r=l " ) 

Center 

55 

Yes 

Yes 

Ball 

(r=l " ) 

Center 

60 

Yes 

Yes 

Ball 

(r=l " ) 

Center 

65 

Yes 

Yes 

Ball 

(r=l " ) 

Center 

70 

Yes 

Yes 

Ball 

(r=l " ) 

On  stiffener 

10 

Yes 

Yes 

Ball 

(r=l " ) 

On  stiffener 

30 

Yes 

Yes 

Ball 

(r=l " ) 

On  stiffener 

35 

Yes 

Yes 

Ball 

(r=l " ) 

On  stiffener 

55 

Yes 

Yes 

Ball 

(r=l " ) 

On  stiffener 

60 

Yes 

Yes 

Ball 

(r=l " ) 

On  stiffener 

65 

Yes 

Yes 

Ball 

(r=l " ) 

On  stiffener 

70 

Yes 

Yes 
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IV.  RESULTS/DISCUSSIONS 


A.  DYNAMIC  IMPACT  SIMULATION 

The  impactor  chosen  for  dynamic  impact  simulation  was  a 
0.0254  cm  (1  in.)  radius  rigid  ball.  The  rigid  ball 
impacted  a  shell  panel  at  various  velocities  and  locations 
(Table  3)  creating  different  dent  sizes  and  depths.  As 
expected,  the  dent  sizes  and  depths  increased  with  the 
increasing  impact  velocity.  Dent  sizes  varied  from  6.35  cm 
(2.5  in.)  to  slightly  less  than  15.24  cm  (6.0  in.)  in 
diameter  as  shown  in  Table  4 .  The  depth  of  dents  varied 
from  about  a  0.3175  cm  (1/8  in.)  to  a  little  over  6.35  cm 
(2.5  in.) .  Although  dent  sizes  and  depths  were  not  the  same 
for  the  deformed  stiffened  and  unstiffened  panels  subjected 
to  the  same  impact  condition,  both  did  show  the  same  trend 
of  an  increase  in  dent  sizes  and  depths  with  respect  to  an 
increase  in  impact  velocity.  At  lower  impact  velocities, 
however,  it  was  difficult  to  determine  the  dent  radius 
because  the  dent  sizes  were  too  small  to  measure  accurately 
with  finite  element  mesh  size  of  0.635  cm  (0.25  in)  .  For 
instance,  with  stiffened  panel  impacted  at  the  velocity  of 
10  m/s  at  the  center,  the  dent  size  was  probably  less  than 
0.635  cm  and  therefore  too  small  to  be  measured.  The 
existence  of  a  dent  was  confirmed  by  the  contour  plot  of 
equivalent  plastic  strain. 
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Table  4  . 


Fuselage  Dent  Results 


Fuselage  Simulation  Dent  Result  Summary 


Panel  Type 

Impact 

Location 

Impact 

Speed 

(m/s) 

Dent 

Radius 

(cm) 

Dent 

Radius 

(in) 

Dent 

Depth 

(m) 

Dent 

Depth 

(in) 

CJndeformed  Panel 

N/A 

N/A 

Undeformed  Stiffened  Panel 

N/A 

N/A 

Deformed  Unstiffened  Panel 

Center 

10 

3.18 

1.25 

0.0030 

0.12 

Center 

30 

6.67 

2 . 98 

0.0152 

0.60 

Center 

35 

7.58 

2 . 60 

0.0200 

0.79 

Center 

55 

6.61 

2.58 

0.0314 

1.24 

Center 

60 

6.55 

2.76 

0.0375 

1.48 

Center 

65 

7.00 

2 . 66 

0.0418 

1 . 65 

Center 

70 

6.75 

0.00 

0.0469 

1 . 85 

Deformed  Stiffened  Panel 

Center 

10 

0.0027 

0.11 

Center 

30 

5.06 

2.39 

0.0115 

0.45 

Center 

35 

6.07 

2.70 

0.0141 

0.56 

Center 

55 

6.87 

2 . 69 

0.0189 

0.74 

Center 

60 

6.83 

2.75 

0.0223 

0.88 

Center 

65 

7.00 

2 . 81 

0.0269 

1.06 

Center 

70 

7 . 15 

0.00 

0.0317 

1.25 

On  stiffener 

30 

0.0041 

0.16 

On  stiffener 

35 

3.94 

2.48 

0.0045 

0.18 

On  stiffener 

55 

6.31 

2.23 

0.0248 

0.97 

On  stiffener 

60 

5.66 

2.25 

0.0286 

1 . 12 

On  stiffener 

65 

5.72 

1 . 66 

0.0316 

1.24 

On  stiffener 

70 

4.21 

0.00 

0.0347 

1.36 

Figure  5  shows  the  plot  of  dent  radius  versus  impact 
velocity.  The  plot  showed  that  initially  there  was  an 
increase  in  the  dent  size  with  an  increasing  impact  velocity 
for  all  three  cases.  For  the  center  impact  of  stiffened  and 
unstiffened  panels,  the  dent  sizes  showed  a  large  increase 
from  a  low  velocity  (10  m/s)  to  a  higher  velocity  (30  m/s) 
and  slowed  down  or  stabilized  at  a  much  higher  velocity.  The 
dent  size  of  the  stiffened  panel  with  an  impact  at  10  m/s 
was  too  small  to  measure  even  from  the  contour  plot  of 
equivalent  plastic  strain. 
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For  the  direct  stiffener  site  impact,  however,  the 
figure  showed  the  dent  radius  experienced  a  large  jump 
initially  and  decreased  at  high  velocities.  There  is  no 
clear  explanation  for  that  at  this  time.  This  phenomenon 
needs  to  be  further  investigated. 


60 

50 

40 

30 

20 

10 

0 


Dent  Radius  vs  Impact  Speed 


10  20  30  40  50 

Impact  Speed  (m/s) 


60 


70 


Figure  5 . 


Dent  size  impact  results. 
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Dent  Depth  vs  Impact  Speed 


— Unstiffened  Panel 


-fl — Stiff  Panel 
Center  Impact 

— Stiff  Panel 

Stiffener  Impact 


Figure  6.  Dent  depths  impact  results. 

The  dent  depths  were  determined  by  the  maximum 
transverse  displacement  shown  on  the  displacement  contour 
plot  of  each  deformed  panel.  Figure  6  showed  that  the  dent 
depth  increases  monotonically  as  a  function  of  impact 
velocity.  Both  the  unstiffened  panel  and  the  direct 
stiffener  impact  panel  showed  a  nearly  linear  relationship 
between  the  dent  depth  and  the  impact  speed.  However,  the 
stiffened  panel  with  the  center  impact  indicated  a  less 
linear  relationship. 

B.  EIGENVALUE  ANALYSES  OF  VIRGIN  PANELS 

Eigenvalue  analyses  were  performed  for  all  virgin 
panels  before  impact  in  order  to  compute  their  initial 
buckling  loads.  The  first  five  buckling  modes  and  their 
associated  buckling  loads  were  calculated.  However,  the 
lowest  buckling  load  was  selected  as  the  reference  value  to 
be  compared  with  the  compressive  failure  load  after  the 
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impact  damage.  The  critical  buckling  loads  for 
unstiffened  and  stiffened  virgin  panels  were  990  N 
15, 920  N,  respectively.  Some  of  the  buckling  mode  shapes 
unstiffened  panels  are  shown  in  (Figure  7-10)  . 


the 

and 

for 
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Step:  Step-1,  Buckling-  prediction  of  undeformed  plate. 
Mode  1:  Eigenvalue  =  989.65 

Primary  Var:  U,  Magnitude 

Deformed  Var:  U  Deformat  ion  Scale  Factor:  +5.080e-02 


Figure  7 . 


Virgin  unstiffened  panel  buckling  Mode 
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Magnitude 
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ODB:  Buckle  Plate.  odk> 


ABAQUS/ STANDARD  Version  6.6-1 


Wed  Mar  28  16:13:01 


Step:  Step-1,  Buckling  prediction  of  unde formed  plate. 
Mode  2:  Eigenvalue  =  2759.4 

Primary  Var:  U,  Magnitude 

Deformed  Var:  U  Deformation  Scale  Factor:  +5.080e-02 


Figure  8.  Virgin  unstiffened  panel  buckling  Mode  2. 
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U,  Magnitude 
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ODB:  Buckle_Plat e . odto  ABAQUS/ STANDARD  Version  6.6-1  Wed  Mar  28  16:13:01  : 

Step:  Step-1,  Buckling  prediction  of  undeformed  plate. 

Mode  3 :  Eigenvalue  =  4069 . 5 

Primary  Var:  U,  Magnitude 


Deformed  Var:  U  Deformation  Scale  Factor:  +5.080e-02 


Figure  9. 


Virgin  unstiffened  panel  buckling  Mode  3. 


U,  Magnitude 
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ODB: 


Buckle_Plate. odto  ABAQUS/ STANDARD  Version  6.6-1  Wed  Mar  28  16:13:01  : 


Step:  Step-1,  Buckling  prediction  of  undeformed  plate. 
Mode  4:  Eigenvalue  =  5908.5 

Primary  Var:  U,  Magnitude 

Deformed  Var:  U  Deformation  Scale  Factor:  +5.080e-02 


Figure  10.  Virgin  unstiffened  panel  buckling  Mode  4. 
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C.  COMPRESSION  TEST/POSTBUCKLING  SIMULATION 

Following  the  dynamic  impact  analyses,  the  deformed 
shell  panels  including  the  residual  stresses  and  strains 
were  imported  into  a  new  model  and  were  subjected  to  in¬ 
plane  compressive  loading.  The  compressive  failure  loads 
from  the  compression  analyses  were  determined  from  the  load 
and  displacement  data.  For  the  unstiffened  panels,  critical 
failure  loads  can  be  read  directly  from  their  load  versus 
displacement  curves.  For  instance.  Figure  11  showed  load 
versus  displacement  curve  for  a  deformed  plate  formed  by  a 
35  m/s  impact.  The  curve  clearly  showed  the  critical  load 
at  750  N.  The  failure  modes  for  unstiffened  panels  were  as 
expected  (Figure  12) .  Their  shapes  were  a  half  sine  curve. 

The  failure  loads  for  stiffened  panels,  on  the  other 
hand,  were  not  as  obvious  by  looking  at  the  load 
displacement  curve  (Figure  13)  .  In  this  case,  the  failure 
loads  were  determined  from  the  last  data  point  in  the  graph 
where  the  finite  element  solution  diverges. 
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Load  vs  Displacement  (N41) 

(Deformed  Plate  @V=35m.'s) 
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Figure  1 1 . 


Load  displacement  for  deformed  plate  at 
V=35m/ s . 


Figure  12 . 


Shell  panel  after  impact  (35m/s)  and 
compression  tests. 
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Load  vs  Displacement  (N  536) 

(Stiffened  Plate  Center  Impart®  V=55m/s) 


In-Plane  Disp  Edge  Node  (m) 


Figure  13.  Load  displacement  for  deformed  stiffened 

plate  at  V=55  m/s. 


Figure  14.  Shell  panel  after  impact  (55m/s)  and 

compression  tests. 


The  compressive  test  analyses  for  the  damaged 
unstiffened  panels  showed  that  certain  fuselage  dents  seemed 
to  increase  the  failure  load  of  the  panel  compared  to  the 
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buckling  load  of  the  virgin  panel.  An  experimental  study  by 
Guijt  also  noted  a  strengthening  effect  caused  by  existence 
of  dents  (Ref.  1)  .  From  the  summary  listed  in  Table  4,  it 
can  be  seen  that  at  low  impact  velocities,  the  compressive 
failure  loads  were  smaller  than  the  virgin  buckling  load. 
In  fact,  at  the  impact  speed  of  lOm/s,  30m/s,  and  35m/s,  the 
compressive  failure  loads  were  606  N,  659  N,  and  725  N, 
respectively,  which  were  less  than  that  of  the  virgin  panel 
(990  N)  .  As  the  impact  velocity  was  increased,  the  failure 
load  continued  to  increase  exceeding  the  buckling  load  of 
the  virgin  panel.  The  failure  load  increased  to  a  maximum 
at  1368  N  and  began  to  decrease.  A  plot  of  failure  load 
versus  impact  speed  is  given  in  Figure  13.  The  failure 
loads  have  been  normalized  against  the  buckling  load  of  the 
virgin  panel.  The  plot  showed  a  general  stiffening  effect 
caused  by  impacts  at  higher  velocities.  An  explanation  of 
this  stiffening  effect  is  given  in  the  last  section  of  this 
chapter . 

Table  5.  Unstiffened  Panel  Buckling  Load. 


Impact 

Location 

Impact 

Speed 

(m/s) 

Dent 

Radius 

(cm) 

Dent 

Radius 

(in) 

Dent 

Depth 

(cm) 

Dent 

Depth 

(in) 

Failure 
Load  (N) 

Center 

10 

3.18 

1.25 

0.30 

0.12 

606 

Center 

30 

6.67 

2 . 63 

1.52 

0.60 

659 

Center 

35 

7.58 

2 . 98 

2.00 

0.79 

725 

Center 

55 

6.61 

2 . 60 

3.14 

1.24 

1093 

Center 

60 

6.55 

2.58 

3.75 

1.48 

1206 

Center 

65 

7.00 

2.76 

4 .18 

1 . 65 

1368 

Center 

70 

6.75 

2 . 66 

4 . 69 

1 . 85 

1351 
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Buckling  Load  at  Various  Impact  Speed 
(Unstiffened  Panel) 


Figure  15.  Buckling  load  for  deformed  unstiffened  panel. 

For  the  deformed  stiffened  shell  panels,  results  from 
compression  test  analyses  revealed  similar  behavior  as 
earlier.  The  failure  loads  for  damaged  stiffened  panels 
started  out  at  small  values  and  increased,  surpassing  the 
buckling  load  for  the  virgin  panel  (15,920  N) ,  to  a  maximum 
and  began  to  decreased  back  down.  From  Table  7,  for  the 

stiffened  panel  with  center  impact  at  10  m/s,  the  failure 

load  was  15,190  N.  This  was  slightly  lower  than  buckling 
load  for  the  virgin  panel.  The  failure  load  reached  a 
maximum  value  of  42,260  N  at  30  m/s  and  started  to  decrease 
rapidly  to  15967  N  at  70  m/s .  These  results,  again, 

appeared  to  suggest  a  stiffening  effect  caused  by  high 
velocity  impacts. 

For  stiffened  panels  with  direct  damage  to  the 

stiffener  site,  it  was  found  that  the  maximum  buckling  load 
was  39,  230  N  at  the  impact  velocity  of  55  m/s .  At  low 
impact  velocities,  the  buckling  loads  were  approximately  25 

%  lower  than  that  of  the  virgin  panel  without  dent. 
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It 


appeared  that  damage  to  the  stiffener  support  at  low 
velocity  impact  can  be  quite  detrimental  to  the  structure  of 
a  fuselage.  The  failure  load  versus  impact  speed  plot  in 
Figure  16  and  17  below  showed  a  summary  of  buckling  behavior 
of  deformed  stiffened  panels. 


Table  6.  Stiffened  Panel  Buckling  Load  Results 

Summary . 


Impact 

Location 

Impact 

Speed 

(m/s) 

Dent 

Radius 

(cm) 

Dent 

Radius 

(in) 

Dent 

Depth 

(cm) 

Dent 

Depth 

(in) 

Buckling 
Load  (N) 

Center 

10 

0.27 

0.11 

15186 

Center 

30 

5.06 

2.39 

1 . 15 

0.45 

42261 

Center 

35 

6.07 

2.70 

1.41 

0.56 

39295 

Center 

55 

6.87 

2 . 69 

1 .89 

0.74 

18044 

Center 

60 

6.83 

2.75 

2.23 

0.88 

16278 

Center 

65 

7.00 

2 . 81 

2 . 69 

1.06 

17690 

Center 

70 

7 . 15 

0.00 

3.17 

1.25 

15967 

On  stiffener 

30 

0.41 

0.16 

11459 

On  stiffener 

35 

3.94 

2.48 

0.45 

0.18 

12084 

On  stiffener 

55 

6.31 

2.23 

2.48 

0.97 

39232 

On  stiffener 

60 

5.66 

2.25 

2 . 86 

1 . 12 

32719 

On  stiffener 

65 

5.72 

1.66 

3.16 

1.24 

28121 

On  stiffener 

70 

4.21 

0.00 

3.47 

1.36 

28719 

Figure  16. 


Buckling  loads  for  stiffened  panel  center 
impact . 


Buckling  Load  at  Various  Impact  Speeds 
(Stiffened  Panel  Stiffener  Impact) 


Figure  17.  Buckling  loads  for  stiffened  panel  stiffener 

impact . 

D.  EIGENVALUE  ANALYSES  OF  DAMAGED  PANELS 

Linear  eigenvalue  analyses  were  performed  for  each 
damaged  fuselage  panel  after  impact  in  order  to  compute 
their  buckling  loads.  The  geometric  models  of  the  damage 
panels  were  imported  from  the  impact  analysis  models  for 
eigenvalue  analysis.  Because  the  linear  eigenvalue  analysis 
does  not  require  residual  stress/strain,  this  information 
was  not  imported.  The  buckling  loads  of  the  damaged  panels 
were  compared  with  the  compressive  failure  loads  after  the 
impact  damage.  Even  though  the  boundary  conditions  and  the 
loading  direction  were  the  same  for  the  two  analyses  (i.e. 
nonlinear  compressive  and  linear  eigenvalue  analyses) ,  the 
comparisons  showed  a  large  difference  in  the  failure  loads 
and  mode  shapes  between  the  two  analyses  as  seen  in  Table  7 . 
For  instance,  at  70  m/s  impact  on  stiffener,  the  damaged 
stiffened  panel  showed  a  buckling  load  of  15600  N  while  the 
failure  load  from  compressive  test  analysis  gave  slightly 
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over  28700  N.  Moreover,  when  comparing  buckling  mode  shapes 
for  the  same  damaged  panel,  the  compressive  test  analysis 
showed  a  larger  deformation  along  the  stiffener  support 
(Figure  18)  .  On  the  other  hand,  the  linear  eigenvalue 
analysis  showed  a  larger  deformation  on  one  of  the  free 
edges  (Figure  19)  .  The  next  four  eigenvalue  mode  shapes 
for  the  same  panel  were  also  different  from  the  compressive 
test  analysis.  The  results  for  the  compressive  failure 
analysis  and  for  the  first  mode  eigenvalue  calculation  are 
shown  in  Figure  18  and  19  respectively. 


S,  Mises 

SNEG,  ( f ract ion  =  -1.0) 
(Avg:  75%) 

+5.687e+08 
+5 . 222e+08 
+4 . 757e+08 
+4.292e+08 
+3. 826e+08 
+3. 361e+08 
+2. 896e+08 
+2 . 4  31e+0  8 
+1 . 965e+08 
+1. 500e+08 
+1 . 0  35e+08 
+5.695e+07 
+1 . 042e+07 


Postbuckling  analysis  of 
ODB:  FOSTBUCKLE_Stif 

Step:  Step-1 

■  Increment  100:  Arc  Length  =  13. 

Primary  Var:  S,  Mises 
Deformed  Var:  U  Deformation  Scale  Factor 


speed 

Sun  Mar  04  13: 


Figure  18.  Buckled  stiffened  panel  with  stiffener  impact 

at  70m/s  (Compressive  Test  Analysis) . 
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U,  Magnitude 
+1 . 000e+00 
+9. 171e-01 
+8. 337e-01 
+7. 504e-01 
+6.670e-01 
+5. 836e-01 
+5.002e-01 
+4. 169e-01 
+3. 335e-01 
+2. 501e-01 
+1.667e-01 
+8. 337e-02 
+0 . 000e+00 


ODB: 


Step:  Step-1 

1  Mode  1:  Eigenvalue  =  15611. 

Primary  Var:  U,  Magnitude 


inpact  V=70m/ s ;  e 
ion  6.6-3  Sun  Mar  04  14:09:26 


Deformed  Var:  U  Deformation  Scale  Factor:  +5.113e-02 


E 


Figure  19. 

Table  7  . 


Results  from  Eigenvalue  Calculation  Mode  1 . 
Compare  Eigenvalues  and  Buckling  Loads. 


Fuselage  Simulation  Eigenvalues  Result  Summary 

Panel  Type 

Impact  Type 

Impact 

Location 

Impact 

Speed 

(m/s) 

Mode  1 
Eigenvalue 

Buckling 
Load (N) 

CJndeformed  Panel 

N/A 

N/A 

N/A 

990 

Undeformed  Stiffened  Panel 

N/A 

N/A 

N/A 

15922 

Deformed  Unstiffened  Panel 

Ball  ( r=l " ) 

Center 

10 

2125 

606 

Ball  ( r=l " ) 

Center 

65 

7513 

1368 

Deformed  Stiffened  Panel 

Ball  ( r=l " ) 

Center 

10 

8978 

15186 

Ball  ( r=l " ) 

Center 

25 

20094 

13770 

Ball  ( r=l " ) 

Center 

55 

20687 

18044 

Ball  ( r=l " ) 

Center 

60 

23357 

16278 

Ball  ( r=l " ) 

Center 

65 

18858 

17690 

Ball  ( r=l " ) 

Center 

70 

6245 

15967 

Ball  ( r=l " ) 

On  stiffener 

70 

15611 

28719 

E.  EXPLANATION  FOR  STIFFENING  EFFECT  OF  FUSELAGE  DENTS 

From  the  compressive  test  analyses,  it  was  noticed  that 
the  existence  of  a  dent  seemed  to  cause  a  strengthening  of 
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the  deformed  fuselage  panels  as  discussed  in  previous 
sections.  It  was  noticed  that  at  lower  impact  velocity,  the 
failure  load  was  smaller  than  that  of  the  undented  virgin 
panel.  At  a  higher  velocity,  the  failure  load  rises  above 
that  of  the  virgin  panel  and  reached  a  maximum  at  a  critical 
impact  velocity.  At  an  impact  velocity  greater  than  the 
critical  velocity,  the  buckling  load  began  to  decrease. 
This  important  observation  was  made  in  both  the  unstiffened 
and  stiffened  panels. 

By  examining  the  dent  sizes  and  depths,  the  dent  shape 
was  not  significant  at  a  low  impact  velocity.  The  dent  was 
a  very  small  localized  deformation  at  the  impact  site 
compared  to  the  rest  of  the  panel.  However,  the  local 
residual  stress  at  the  dent  caused  by  the  impact  reduced  the 
overall  strength  of  the  panel.  As  a  result,  a  very  small 
local  dent  reduced  the  compressive  failure  load  lower  than 
the  buckling  load  of  virgin  plate  without  a  dent.  As  the 
impact  velocity  increased,  the  dent  shape  became  more 
significant.  Here,  there  were  two  competing  factors.  First 
the  dented  site  has  a  lower  local  strength  with  high 
residual  stresses  which  would  reduce  the  compressive  failure 
load.  On  the  other  hand,  the  shape  of  the  dent  played  a 
local  bending  stiffening  effect  which  would  increase  the 
compressive  failure  strength.  A  reduction  in  compressive 
failure  load  can  be  viewed  as  a  negative  effect  while  an 
increase  in  compressive  failure  loads  a  positive  effect. 
Thus,  up  to  a  certain  dent  size  and  depth,  the  net  effect 
(i.e.  negative  effect  subtracted  from  the  positive  effect) 
increased  positively.  As  a  result,  the  compressive  failure 


load  of  a  dented  panel  increased  with  the  impact  velocity. 

The  compressive  failure  load  eventually  exceeded  that  of  the 
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virgin  panel.  After  reaching  a  peak  compressive  failure 
load,  the  net  effect  would  start  to  decrease  resulting  in 
decrease  of  the  compressive  failure  load. 

In  an  attempt  to  support  the  above  explanation,  a 
series  of  finite  element  models  were  created.  The  fuselage 
panel  with  protruding  structure  model  was  described  in  the 
Chapter  2,  Sec.  C.  This  model  was  an  attempt  to  investigate 
the  behavior  of  a  dented  fuselage  panel  via  eigenvalue 
analysis.  The  height  of  protruding  structure  modeled  the 
dent  depth  while  the  size  of  the  cutout  modeled  the  dent 
size.  A  series  of  eigenvalue  analysis  were  performed  using 
various  protruding  heights.  Table  8  summarized  the  results. 

In  the  model,  a  cutout  hole  represents  the  reduced 
strength  of  the  dented  section  caused  by  residual  stresses 
while  the  protrusion  represents  the  increase  of  bending 
stiffness  resulting  from  a  dent.  As  seen  in  the  table,  when 
the  effect  of  the  hole  is  greater  than  the  protrusion 
effect,  the  buckling  load  is  lower  than  that  of  the  panel 
without  a  hole.  However,  as  the  protrusion  effect  (i.e. 
increase  of  bending  stiffness  with  dent)  becomes  more 
dominant,  the  buckling  load  becomes  much  greater  than  that 
of  the  plate  without  a  hole. 
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Table  8  . 


Buckling  Loads  for  panel  with  protrusion. 


Panel  Type 

Protrusion 
Height  h 
(cm) 

Protrusion 
Angle  a 

Buckling 
Load  (N) 

Panel  without  cutout 

N/A 

N/A 

990 

8"  cutout  w/o  protrusion 

0 

0 

660 

8"  cutout  w/  protrusion 

0.1 

8.6 

724 

8"  cutout  w/  protrusion 

0.1 

4.3 

777 

8"  cutout  w/  protrusion 

0.2 

16.8 

800 

8"  cutout  w/  protrusion 

0.3 

24 . 4 

874 

8"  cutout  w/  protrusion 

0.635 

90 

967 

8"  cutout  w/  protrusion 

1 . 0 

60.9 

1047 

8"  cutout  w/  protrusion 

1 . 0 

40.8 

1078 

8"  cutout  w/  protrusion 

1 . 0 

24.8 

1145 

8"  cutout  w/  protrusion 

2.0 

59.9 

3272 
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V.  CONCLUSION/RECOMMENDATION 


This  thesis  investigated  the  effect  of  fuselage  dents 
on  compressive  failure  load  using  computer  modeling  and 
simulations.  By  modeling  the  different  impact  velocities 
and  locations,  various  dent  shapes  (sizes  and  depths)  were 
created  for  both  unstiffened  and  stiffened  panels.  An 

additional  study  was  also  performed  to  attempt  to  further 
explain  the  stiffening  effect  of  fuselage  dents  on  the 

compressive  failure  load.  From  this  study,  a  few 

conclusions  and  generalizations  were  made  based  on 

examination  of  the  results. 

1 .  Depending  on  the  dent  status  in  a  panel  caused  by 
an  impact,  the  dent  may  decrease  or  increase  the 
failure  load  of  the  panel  compared  to  the  virgin 
panel  without  impact. 

2.  It  was  observed  that  at  a  low  impact  velocity, 

failure  loads  of  damaged  panels  were  generally 
smaller  than  that  of  the  virgin  panels.  However, 
as  the  impact  velocity  was  increased,  the  failure 
load  also  increased  and  eventually  surpassed  the 
buckling  load  of  the  virgin  panels.  At  a  certain 
critical  impact  velocity,  the  failure  load  reached 
a  maximum  after  which  it  began  to  decrease.  In 
general,  the  existence  of  a  dent  can  strengthen  a 
deformed  panel  in  a  certain  dent  size  range  or 
impact  velocity. 

3.  A  direct  impact  or  damage  to  the  stiffener  site  at 
a  low  impact  velocity  can  be  detrimental  as  the 
failure  load  was  shown  to  be  much  less  than  that 
of  the  virgin  panel  without  impact. 

4.  The  linear  eigenvalue  calculations  of  the  damaged 
panels  did  not  compared  well  to  the  compressive 
failure  load  of  the  same  panel  because  residual 
stresses  were  not  accounted  in  the  eigenvalue 
analysis . 
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A  few  recommendations  can  be  made  for  future  studies 
based  on  the  results  obtained  in  this  investigation: 

1.  In  reality,  dents  can  be  of  many  geometric  shapes 
and  sizes.  Recommend  using  different  impact 
shapes  (i.e.  a  cylinder  impactor  at  an  angle) . 

2  Only  two  impact  locations  were  investigated  in 
this  study.  Recommend  investigating  other  impact 
locations  to  fully  understand  the  effect  of  dents 
on  stiffened  panel. 

3.  A  hat-like  stiffener  model  was  used  in  this  study. 
Based  on  literature  reviews,  other  stiffener 
designs  are  available  for  fuselage  panel. 
Recommend  using  a  different  stiffener  design. 
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APPENDIX  A. .  STIFFENED  AND  UNSTIFFENED  PANEL  IMPACT 

MODELING  USING  ABAQUS/CAE 


This  appendix  explains  in  detail  how  impact  modeling 
was  created  using  ABAQUS/CAE.  The  final  result  would  be  an 
input  file  which  is  used  to  submit  into  a  job  request  for 
analysis  using  ABAQUS/Explicit .  ABAQUS  is  an  advanced 
finite  element  analysis  that  provides  complete  and  powerful 
solutions  for  linear  and  nonlinear  engineering  problems.  It 
is  a  suite  of  finite  element  analysis  modules.  ABAQUS/CAE, 
having  a  modern  graphical  user  interface  (GUI)  of  menus, 
icons,  and  dialog  boxes,  provides  the  most  complete 
interface  with  the  ABAQUS  solver  programs  available  [Ref  2] . 
In  order  to  understand  the  process,  the  user  should  be 
familiar  with  the  components  of  the  ABAQUS/CAE  and  the 
appearance  of  the  window.  Figure  20  shows  the  components 
that  appear  in  the  main  window. 
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Tillebar  Menubar  Toolbar  Context  bar 


Figure  20.  Components  of  the  ABAQUS/CAE  main  window. 


1.  The  part  (stiffened  plate)  for  the  model  can  be 
created  in  the  part  module  of  the  main  window. 
From  the  main  menu  bar: 

Part  -4  Create  -4  Name  the  part  "Stif fened_Plate" 
■4  Choose  3-D  modeling  space,  deformable  type, 
shell  shape  feature,  and  extrusion  type  feature  -4 
Set  approximate  size  to  1.0  4  Click  OK  (Figure 
21) 
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Create  Part 


Name:  |Stiffened_Plate 

Modeling  Space 

(•  3D  r  2D  Planar  C  Axisymmetric 


Type 

(•  jDeformabje: 
Discrete  rigid 
Analytical  rigid 

Base  Feature 
Shape  — |  r  Type 

Solid 
<?  Shell 
<~  Wire 
C  Point 


Options 


None  available 


|  Planar 

Revolution 

Sweep 


Approximate  siz 

a:  |l 

Continue. . . 

Cancel 

Figure  2 1 . 


Create  Part  Dialog  Box 


2.  From  the  edit  tools  to  the  left  of  the  viewport, 
select  the  "+"  (create  isolated  point  tool)  to 
enter  following  coordinates:  (-0.254,0)  and 

(0.254,0)  .  Select  "create  lines  tool"  (to  the 
left  of  create  isolated  point  tool)  and  connect 
the  two  points  previously  created.  Follow  the 
specification  given  in  Table  2  to  create  two  hat 
stiffeners  at  0.2032  m  (8  inches)  apart  using  the 
"create  isolated  point"  and  "create  lines"  tools. 
Once  complete  (see  Figure  22) : 

Click  Done  at  the  bottom  of  viewport  ->  in  Edit 
Base  Extrusion  box,  enter  0.508  for  depth  4  Click 
OK  (Figure  23) 
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Line  part  shell  extrusion 


Unmeshed  Stiffened  Panel 


3.  The  rigid  ball  impact  is  created  in  a  similar 
fashion  as  earlier.  From  the  main  menu  bar: 

Part  ->  Create  3D  Modeling  Space  ->  Analytic 
Rigid  Type  Revolved  shell  Approximate  size  1 
->  Click  Continue .  ->  Select  Create  Arc  tool 

Follow  direction  at  bottom  of  viewport;  select 
center  point  at  origin;  select  perimeter  point  at 
(0,-0.0254)  and  second  point  at  (0.0254,0).  4 

Create  second  arc  with  at  origin  and  perimeter 
points  at  (0.0254,0)  and  (0,0.0254).  ->  Click  OK. 

Create  reference  point  at  top  of  ball  ->  Tools  -> 
Reference  Point  ->  Select  top  point  of  ball  -> 
Click  Done.  See  Figure  24. 


Figure  24.  Rigid  Impact  Ball  with  reference  point. 


37 


The  material  property  is  assigned  to  the  stiffened 
plate  in  the  Property  Module.  First,  the  material 
property  must  be  defined.  From  the  main  menu  bar: 


Material  Create  Name  material  "A1-T3"  In 
General  pull  down  menu,  select  density,  enter  2780 
kg/m3  for  mass  density.  In  Mechanical  pull  down 
menu,  select  Elasticity  ->  Elastic  ->  Enter 
Young' s  Modulus  and  Poisson' s  Ratio  provided  in 
Table  1.  ->  Mechanical  ->  Plasticity  ->  Plastic  -> 
Enter  Yield  Stress  and  Plastic  Strain  in 
accordance  with  Table  1.  Click  OK. 


Edit  Material 


Name:  AI-T3 

Material  Behaviors 


Density 

Elastic 


General  Mechanical  Ihermal  Other 

Plastic 

Hardening:  | Isotropic  ~^\ 

Use  strain-rate-dependent  data 
l~”  Use  temperature-dependent  data 
Number  of  field  variables:  |  0  ~^~j 

Data 


Yield 

Plastic 

Stress 

Strain 

1 

3.47E+008 

0 

2 

5.699E+008 

0.1577 

Delete 


▼  Suboptions 


Cancel 


Figure  25.  Edit  Material  Property  Dialog  Box. 
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Section  -4  Create  4  Select  Shell  Category  -4 
Homogeneous  Type  -4  Click  OK.  -4  In  Edit  Section 
window,  enter  shell  thickness  0.0015875  m  -4 
Select  Material  A1-T3  -4  Click  OK.  (Create  second 
section  "Section-2"  with  shell  thickness  0.003175 
m)  . 

Assign  -4  Section  -4  Select  both  stiffener  joints 
in  viewport  -4  Click  OK  -4  Select  "Section-2". 
Assign  "Section-1"  to  the  rest  of  plate  using  the 
same  procedure . 

5.  Parts  assembly  is  accomplished  next.  This  is 
performed  in  the  Assembly  Module  by  creating  an 
instance  of  each  part  in  model  and  reorient  to  the 
desired  position.  From  main  menu  bar: 

Instance  -4  Create  -4  Select  Rig_Ball  -4  Click 
Apply  -4  Select  Stif fened_Panel 


Figure  26.  Ball  &  Stiffened  Plate  Impact  Analysis  Setup. 
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6.  Analysis  step  is  created  next  via  the  Step  Module. 
This  part  allows  selection  of  type  of  analysis 
needed  to  perform.  In  this  case,  dynamic  explicit 
analysis  is  chosen.  From  the  main  menu  bar: 

Step  Create  ->  General  procedure  type  ->  Dynamic 
Explicit  ->  Continue  ->  In  Edit  Step  window.  Click 
OK. 

7.  Output  can  be  requested  in  the  Step  Module.  For 
this  analysis,  the  default  output  request  is  used. 
To  rename  the  both  the  history  output  and  field 
output  request: 

Output  ->  Field  Output  Requests  or  History  Output 
Requests-^  Rename  Enter  name  or  use  default. 

8.  Contact  and  interaction  property  are  defined  in 
the  Interaction  Module.  From  main  menu  bar: 

Interaction  ->  Property  ->  Create  ->  Contact  type 
->  Continue  ->  Mechanical  -^Tangential  Behavior  -> 
Frictionless  Mechanical  Normal  Behavior  -> 
"Hard"  Contact  Allow  separation  after  impact 
Click  OK. 

Interaction  ->  Create  ->  Surface-to-surface 
contact  (Explicit)  ->  Continue  ->  Follow  direction 
at  bottom  of  viewport  ->  select  plate  as  first 
surface  ->  select  rigid  ball  as  second  surface  -> 
Click  Done . 

9.  Boundary  conditions  and  initial  condition  are 
defined  in  the  Load  Module.  The  plate  is  rigidly 
constrained  on  four  edges.  The  reference  point  on 
the  ball  is  also  rigidly  constrained  except  for 
the  direction  of  motion  of  impact  (U2)  .  The 
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initial  condition  or  velocity  of  Rigid  ball  is 
applied  at  the  ball's  reference  point.  From  the 
Main  Menu  Bar:  (Figure  27) 

BC  Name  Plate_Edge_BC  Select  Step-1  -> 

Category  Mechanical  ->  Types  Displacement/Rotation 
->  Continue  ->  Follow  instruction  at  bottom  of 
viewport  ->  select  all  four  edges  In  Edit 

Boundary  Condition  Box,  select  all  6  DOFs  (3 
displacements/  3  rotations)  ->  Click  OK 

Predefined  Field  ->  Create  ->  Initial  Step  -> 
Category  Mechanical  ->  Types  Velocity  Continue 
->  Select  reference  point  (RP)  on  Rigid  Ball 
Enter  velocity  for  U2  (10,  30,  35,  55,  60,  65,  or 
70m/s) . 


2 


Figure  27.  Applied  Boundary  Condition  &  Initial  Velocity 


10.  In  order  to  mesh  the  fuselage  panel,  it  must  first 
be  seeded.  There  are  various  ways  to  seed,  but 
seed  by  "Edge  Size"  will  be  used  here.  Once  the 
panel  is  seeded,  the  model  can  be  selected  to 


41 


mesh.  Depending  on  what  kind  of  instance  was 
created  in  the  Assembly  Module  (Dependent  or 
Independent),  the  model  must  be  meshed  as  such. 
If  the  instance  was  Dependent,  then  the  model  must 
be  meshed  via  the  original  part  because  the 
instance  depends  on  original  part.  If  instance 
was  Independent,  then  model  must  be  meshed  via 
assembled  part  because  it  does  not  depend  on  the 
original  part.  A  part  modeled  as  rigid  (rigid 
ball)  can  not  meshed  in  ABAQUS.  In  this  case,  we 
seed  the  flat  panel  separately  from  the  stiffener 
support.  In  the  main  menu  bar: 

Seed  Edge  by  size  ->  select  an  edge  of 

stiffener  to  create  local  seed  Click  Done  -> 
Enter  seed  size  0.00635  m  (Figure  28). 


Figure  28.  Edge  of  panel  seeded  at  0.00635  m. 


Seed  ->  Edge  by  size  ->  Select  edge  of  stiffener 
support  to  create  local  ->  Click  Done  ->  Enter 
size  of  seed  0.00635  m  (Figure  29) . 
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Figure  29.  Edge  of  panel  and  stiffener  support  seeded  at 

0 . 00635  m. 

Mesh  ->  Instance  ->  Select  the  whole  panel 
instance  to  mesh  ->  Click  Done  (Figure  30) 


Figure  30.  Meshed  stiffened  panel  at  0.00635  m  seed. 

11.  The  last  step  is  to  create  a  job  to  be  submitted 
for  analysis  using  the  Job  Module.  ABAQUS  allows 
different  types  of  job  submission  such  as  full 
analysis,  restart,  recover  (from  a  crash)  or 
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continue  analysis.  It  also  allows  configuration 
of  computer  memory  usage  for  analysis  as  well  as 
option  for  multiple  processors  to  run  in  parallel. 
From  menu  bar: 

Job  ->  Create  ->  Submission  tab,  select  Full 
analysis  Memory  tab,  enter  512  for  both  -> 
Precision  tab  ->  ABAQUS/Explicit  precision  double 
->  Nodal  out  precision  full  ->  Click  OK.  (Figure 
31) 


Figure  31.  Create  a  job  for  analysis. 

12.  At  this  point  an  input  file  should've  been 
created  and  can  be  submitted  to  analysis.  The 
different  impact  analysis  as  listed  in  Table  3  can 
be  performed  for  unstiffened  panel  by  simply 
changing  the  initial  velocity  and  resubmit  the 
j  ob . 
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APPENDIX  B.  STIFFENED  AND  UNSTIFFENED  PANEL 
EIGENVALUE  CALCULATION  AND  COMPRESSION  TEST 
MODELING  USING  ABAQUS/CAE 


This  appendix  explains  in  detail  how  models  for 
eigenvalue  calculation  analysis  and  compression  test 
analysis  were  created  using  ABAQUS/CAE.  The  final  result 
would  be  an  input  file  which  is  used  to  submit  into  a  job 
reguest  for  analysis  using  ABAQUS/Standard .  ABAQUS  is  an 
advanced  finite  element  analysis  that  provides  complete  and 
powerful  solutions  for  linear  and  nonlinear  engineering 
problems.  It  is  a  suite  of  finite  element  analysis  modules. 
ABAQUS/CAE,  having  a  modern  graphical  user  interface  (GUI) 
of  menus,  icons,  and  dialog  boxes,  provides  the  most 
complete  interface  with  the  ABAQUS  solver  programs  available 
[Ref.  2]  .  In  order  to  understand  the  process,  the  user 
should  be  familiar  with  the  components  of  the  ABAQUS/CAE  and 
the  appearance  of  the  window.  Figure  20  shows  the  components 
that  appear  in  the  main  window. 

A.  EIGENVALUE  CALCULATION/BUCKLING  PREDICTION  ANALYSIS 
MODELING  USING  ABAQUS/CAE 

1.  The  part  (deformed  stiffened  plate)  involved  in 
this  analysis  was  imported  from  the  impact 
analysis  created  in  Appendix  A.  into  a  new  model 
using  the  *import  option.  From  the  main  menu  bar 
create  a  new  model  and: 
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File  ->  Import  ->  Part  ->  File  Filter,  select 
Output  Database  (*.odb*)  ->  select  output  file 
IMPACT_Stif  f enedPlate_C_v35  .  odb  ->  Click  OK  -> 
Select  STIFFENED_PANEL-1  ->  Rename  part  name 
"Stif fened_panel"  ->  Select  "Import  deformed 
configuration."  ->  Click  OK.  (Figure  32) 


Figure  32.  Imported  deformed  fuselage  panel  at  35m/s 

impact  speed. 

2 .  Once  the  deformed  part  has  been  imported,  it  has 
to  be  assembled  in  the  Assembly  Module.  Even 
though  only  one  part  is  involved  and  instance  has 
to  be  created.  Ensure  instance  name  in  this  model 
is  same  as  instance  name  in  the  previous  model. 
From  main  menu  bar: 

Instance  ->  Create  ->  in  Create  Instance  box,  name 
instance  "Stiffened  Panel-1  ->  Click  OK. 
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3.  Next,  an  analysis  step  must  be  created  in  the  Step 
Module.  Here,  the  Linear  Perturbation  Buckle 
procedure  is  used.  Output  of  the  analysis  is  also 
requested  in  this  module.  In  the  main  menu  bar: 


Step  4  Create  ->  Select  Linear  Perturbation 
procedure  type  ->  Choose  Buckle  ->  Click  OK  ->  in 
Edit  Step  box,  select  Subspace  Eigensolver  -> 
Request  5  eigenvalues  ->  Click  OK. 

Output  ->  Field  Output  Request  ->  Edit  ->  F-Output- 
1  ->  Click  OK  to  accept  the  default  request. 

4.  Next,  the  imported  geometry  must  be  defined  as 
initial  condition  for  this  new  model.  In 

addition,  boundary  conditions  and  load  must  also 
be  defined  as  well.  These  are  done  in  the  Load 
Module.  From  the  main  menu  bar: 


Predefined  Field  ->  Choose  Step  "Initial"  -> 
Category  "Other"  ->  Type  "Initial  State"  ->  Click 
Continue  Select  the  fuselage  panel  instance  in 
the  viewport  to  be  assigned  initial  state 
Bottom  of  viewport  click  Done  in  Edit 
Predefined  Field,  enter  the  previous  analysis  job 
name  "IMPACT_Stif fenedPlate_C_V35"  ->  Click  OK. 
(Figure  33) 


47 


Figure  33.  Specify  exact  job  name  from  previous  impact 

analysis . 


BC  ->  Create  ->  Name  Simple  Support  ->  Select  Step- 
1  ->  Category  "Mechanical"  ->  Type 
"Displacement/Rotation"  Click  Continue  -> 
Select  the  far  edge  from  screen  to  applied  BC;  may 
need  to  rotate  and  zoom  in  to  select  all  nodes  on 
this  edge->  at  bottom  of  viewport,  click  Done  -> 
in  Edit  Boundary  Condition  box,  set 
U1=U2=U3=UR2=UR3=0  (simple  support)  Click  OK. 
(Figure  xxx) 

BC  ->  Name  "Compressed  Edge"  ->  Step-1  ->  Category 
"Mechanical"  ->  Types  "Displacement/Rotation"  -> 
Continue  ->  Set  BC  for  closest  edge  on  screen 
U 1 =U 2  =UR2  =UR3 = 0  Click  OK.  (Figure  34) 
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Figure  34 . 


Boundary  condition  for  fixed  end  and  loaded 
end . 


Load  ->  Create  ->  Step-l->  Category  "Mechanical"  -> 
Type  "Shell  edge  load"  ->  Continue  ->  Select  same 
edge  as  the  "Compressed  Edge"  BC  ->  Click  Done  at 
bottom  of  viewport  ->  in  Edit  Load  box,  set  1 
Newton  for  "Magnitude"  ->  Click  OK.  Figure  (35  & 
36) 

5.  Last  step  is  to  create  a  job  for  eigenvalue 
calculation  using  the  Job  Module  as  shown  in 
Appendix  A.  For  this  job,  use  default  job 
settings . 

6.  The  model  can  be  submitted  for  analysis  after  the 
job  is  created. 
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Figure  35.  Apply  unit  load  for  eigenvalue  calculation. 


Figure  36.  Boundary  conditions  and  unit  load  applied. 
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B.  COMPRESSIVE  TEST  ANALYSIS  MODELING  USING  ABAQUS/CAE 

1.  Create  a  new  model.  From  the  main  menu  bar: 

Model  ->  Create  ->  Name 

"Stiff enedPlate_C_V35_PostBuckle"  Click  OK. 

2.  Import  the  deformed  geometry  from  impact  analysis 
using  the  same  procedure  as  previous  model. 

3.  Create  a  part  instance  in  the  Assembly  Module 
using  the  same  procedure  as  previous  model. 

4.  Create  analysis  step-1  in  the  Step  Module.  Accept 
the  default  output  request. 

Step  ->  Create  ->  General  procedure  type  ->  Select 
Static,  Riks  method  ->  Continue  ->  Nlgeom  set  to 
"ON"  Click  OK  (Figure  37) 


Figure  37.  Static  Riks  method  used  in  analysis  step. 
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5.  In  the  Load  Module,  create  boundary  conditions  as 
in  the  Eigenvalue  Calculation  analysis.  The 
imported  deformed  geometry  and  material  state  is 
set  as  initial  similar  to  previous  model.  Initial 
load  used  during  Riks  method  is  50000  Newton.  In 
the  main  menu  bar: 

Load  ->  Name  "Compressive  Load"  Step-1  -> 

Category  "Mechanical"  ->  Type  "Shell  edge  load"  -> 
Continue  ->  Select  closest  edge  in  viewport 
Click  Done  in  Edit  Load  box,  enter  50000  Newton 
for  "Magnigtude"  Click  OK.  (Figure  38) 


Figure  38.  Apply  load  for  Riks  buckling  method. 


6.  Create  job  for  analysis  similar  to  previous  model. 
Use  default  setting  for  this  job. 
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7.  Once  the  job  has  been  created.  It  can  be  submitted 
for  a  full  analysis.  All  compressive  failure 
models  can  be  created  in  a  similar  fashion  by 
following  the  above  listed  steps. 


53 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


54 


LIST  OF  REFERENCES 


1.  Guijt,  C.,  "New  Structural  Guidelines  For  Dent 
Allowables  on  Fuselages,"  (unpublished  work). 

2.  Hibbit,  Karlsson  &  Sorensen,  Inc.,  ABAQUS/CAE  User's 
Manual  Version  6.6. 

3.  R.  J.  H.  Wanhill,  "Aircraft  Corrosion  and  Fatigue 
Damage  Assessment,  NLR  Technical  Publication  TP94401  L, 
National  Aerospace  Laboratory,  Amsterdam,  1994. 

4.  Siegel,  M.,  Gunatilake,  P.,  "Remote  Inspection 
Technologies  for  Aircraft  Skin  Inspection, "  IEEE 
Workshop  on  Emergent  Technologies  and  Virtual  System 
for  Instrumentation  and  Measurement,  Niagara  Falls, 
Ontario,  Canada,  May  15-17,  1997. 

5.  http : / / www . matweb . com,  March  27,  2007. 


55 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


56 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center 
Ft.  Belvoir,  Virginia 

2 .  Dudley  Knox  Library 
Naval  Postgraduate  School 
Monterey,  California 

3.  Young  W.  Kwon,  Code  ME/Kw 

Dept,  of  Mechanical  and  Astronautical  Engineering 
Naval  Postgraduate  School 
Monterey,  California 

4.  Nann  C.  Lang 

Monroe,  North  Carolina 


57 


